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1 Introduction 



In the past five years considerable progress has been made towards a QCD- 
based and model independent description of hadrons containing a heavy- 
flavour quark. This progress has been achieved by the use of the heavy 
mass limit for the heavy quark, in which it is replaced by a static source of 
a colour field. This limit of QCD leads to a well-defined field theory, the 
so-called Heavy Quark Effective Theory (HQET). 

The infinite mass limit from QCD has been used for some time for var- 
ious purposes but the main observation, which triggered the enormous 
development of this field, is that in the infinite mass limit two additional 
symmetries appear, which are not present in full QCD ||. Since then the 
implications of the heavy mass limit and HQET have been extensively studied 
in innumerable publications, and the development of the field is documented 
in more or less extensive reviews ||. 

The two additional symmetries of the heavy mass limit have important 
phenomenological applications; they lead to model-independent relations be- 
tween form factors describing e.g. exclusive weak decays. The origin of the 
new symmetries is quite simple. The first symmetry is a heavy-flavour sym- 
metry, which is due to the fact that the interaction of the quarks with the 
gluons is flavour-blind and in the heavy mass limit all heavy quarks act as 
a static source of colour. Formally this corresponds to an SU (2) symmetry 
relating b and c quarks moving with the same velocity. The second symmetry 
is the spin symmetry of the heavy quark. The interaction of the heavy quark 
spin with the "chromomagnetic" field is inversely proportional to the heavy 
mass and hence vanishes in the infinite mass limit. As a consequence, the 
rotations for the heavy quark spin become an SU(2) symmetry, which holds 
for a fixed velocity of the heavy quark. 

Corrections to the limit vhq — + oo may be studied systematically in 
the framework of HQET. The corrections are given as power series expan- 
sions in two small parameters. The first small parameter is the strong cou- 
pling constant, taken at the scale of the heavy quark a s {mo). This type 
of correction may be calculated systematically using perturbation theory in 
HQET. The second type of correction is characterized by the small param- 
eter R/rriQ, where A is a scale of the light QCD degrees of freedom, e.g. 
A ~ rrihadron — itiq- In the effective theory approach this type of corrections 
enter through operators of higher dimension, the matrix elements of which 
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have to be parametrized in general by additional form factors. 

More or less parallel to the development of heavy quark symmetry and 
HQET, which is well suited for exclusive decays, the heavy mass expansion 
has been applied also to inclusive decays [}£|-||I0|. The heavy quark mass sets 
a scale that is large compared to Aq C d, and one may use a similar setup as in 
deep inelastic scattering for the inclusive decays; in particular, the operator 
product expansion for the inclusive decays yields an expansion in inverse 
powers of the heavy quark mass. In this way one may not only study total 
rates, but also differential distributions such as the lepton energy spectra in 
inclusive semileptonic decays. 

This mini-review is intended to set up the stage for some of the contri- 
butions to this conference, which present in some detail new calculations in 
HQET or new results for inclusive heavy hadron decays. It is divided into two 
parts, one dealing with HQET and exclusive decays and the second devoted 
to inclusive decays. 

In the next section we discuss the heavy mass limit of QCD and the 
additional symmetries of this limit. We formulate HQET, including terms 
up to order 1/m 2 . The strategy of a HQET calculation is outlined and applied 
to the weak decay matrix elements relevant for the semileptonic decays B — > 
D^lv. In section 3 we consider the setup for the heavy mass expansion for 
inclusive decays and study the rate and differential distributions in inclusive 
semileptonic decays. Finally we comment on inclusive non-leptonic processes 
and inclusive rare decays and conclude. 



2 Heavy Quark Effective Theory 

The Green functions of QCD containing a heavy quark in general depend on 
its mass niq. This mass sets a scale that is large compared to the scale A, 
which characterizes the light degrees of freedom, A/tuq is small and becomes 
a reasonable expansion parameter. The leading order in this parameter cor- 
responds to the infinite mass limit of QCD, which corresponds to an effective 
theory where the degrees of freedom related to this large scale have been 
removed. This effective theory, the so-called HQET, may be formulated as 
a Lagrangian field theory, and its Lagrangian may be obtained from QCD. 
There are several ways to construct this Lagrangian and the one closest to 



the idea of "integrating out" heavy degrees of freedom is discussed in [O], 
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where the small components of the heavy quark spinor field are identified as 
the heavy degrees of freedom and are removed by integrating over them in 
the generating functional of QCD Green functions. 

We shall not go through any derivation here, but rather state the result 
and its relation to full QCD. We denote the heavy quark field of full QCD 
by Q and define 

Q v (x) = exp(im Q vx)Q(x) = h v (x) + H v {x), (1) 

where v is a velocity (v 2 = 1), which is later identified with the velocity of 
the heavy hadron. Extracting this phase factor from the full QCD field Q 
removes the dominant part vhqv of the heavy quark momentum, since this 
phase redefinition corresponds to a splitting of the heavy quark momentum 
according to p — vhqv + k, where the residual momentum k is small, of 
the order of A. Furthermore, h v (H v ) is the large (small) component field, 
corresponding to the projections 

h v = P + Q v , H V = P.Q V , P± = i(l±#. (2) 

The small component field H v is related to the large scale rriQ] integrating 
out H v from the generating functional of QCD Green functions corresponds 
to the replacement 

H v = P- (n ^-r) 0hv (3) 



,2m + ivD 

and this yields a non-local "Lagrangian" of the form fTT 



C = h v (ivD)h v + h v iI?P-(- -—)il?h v , (4) 

which still contains all orders in I/tuq. However, the non-locality appearing 
in the second term of (Q) may be expanded into an infinite series of local 
terms, which come with increasing powers of 1/mq. Hence one may in this 
way establish the desired heavy mass expansion for the Lagrangian. The first 
few terms of this expansion are 



C = h v (ivD)h v + (—) h v iDPWh v + ( —) h v iDP(-ivD)iDh v + ■ 
\2mJ \2mJ 

(5) 
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The non-local expression (H) is still equivalent to full QCD; in particular it 
is independent of the still arbitrary velocity vector v. In fact, the Lagrangian 
is invariant under an infinitesimal shift of the velocity 

v — > v + 5v v ■ 8v = (6) 

h v ^h v + ^-(l + P- ^-^nA K 

2 V 2m + ivD J 

iD —>■ —m5v. 



This invariance is the so-called reparametrization invariance |12| , which has 
non-trivial consequences for the Lagrangian and also for matrix elements, 
since it relates terms of different orders of the 1/m expansion. 

However, the increasing powers of 1/ttiq have to be compensated by the 
dimension of the operators appearing in the expansion. In a field theory, these 
operators are not a priori defined, since they have to be renormalized. This 
renormalization leads to additional dependences on the heavy mass, which 
are in general logarithmic, at least in perturbation theory. The expansion of 

thus gives only the coefficients of the operators at the scale tuq, at which 
the heavy degrees of freedom are integrated out. 

The leading logarithmic corrections to the Lagrangian have been calcu- 



lated |T3| and the result at some scale < tuq is 
C = C h v (tvD)h v (7) 



+ {—) [CoCMtDYK - C C 2 h v (ivD) 2 h v + (-t)C C 3 h v a^tD»iD u h v 
+ ■■■ 

with the coefficients (in Feynman gauge) 

C = r7 -8/(33-2n / ) ^ = j (g) 
C 2 = 3^-8/(33-2^) _ 2 = v -9/(33-2n f ) 

where 

v - 



a s (m Q )' 

The fact that C% = 1 is a consequence of reparametrization invariance, 
which implies non-trivial relations between the renormalizations of the var- 
ious terms in the Lagrangian; e.g. some of the renormalization constants 
of the second-order terms in the Lagrangian may be calculated from the 



first-order ones JT1] 
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2.1 The Heavy Quark Limit and Additional Symme- 
tries 

The leading term of the Lagrangian defines the heavy-quark limit and 
exhibits the two additional symmetries. The first symmetry is the heavy- 
flavour symmetry, which is due to the fact that the leading term in the 
heavy mass expansion is mass-independent. If there are two heavy flavours 
described by the operators b v and c v , then the total Lagrangian is simply the 
sum of the two 

£ = b v {ivD)b v + c v (ivD)c v , (9) 
which is invariant under SU(2)hf rotations among the two fields 

Uv(fy UeSU(2) HF . (10) 

We have put a subscript v for the transformation matrix U, since this sym- 
metry only relates heavy quarks if they move with the same velocity. In other 
words, there is a heavy-flavour symmetry in each velocity sector. 

The second symmetry is the heavy-quark spin symmetry. As is clear form 
the Lagrangian in the heavy-mass limit, both spin degrees of freedom of the 
heavy quark couple in the same way to the heavy quark; we may rewrite the 
leading-order Lagrangian as 

C = h + v s {ivD)h+ s + K s (ivD)h- s , (11) 

where hf s are the projections of the heavy quark field on a definite spin 
direction s 

h^ s = l -{l± 1 4)h v , s-v = 0. (12) 

This Lagrangian has a symmetry under the rotations of the heavy quark spin 
and hence all the heavy hadron states moving with the velocity v fall into 
spin-symmetry doublets as itlq — > oo. In Hilbert space this symmetry is 
generated by operators S v (e) as 

[h v , S v (e)] = i^ 5 h v (13) 

where e with e 2 = —1 is the rotation axis. The simplest spin-symmetry 
doublet in the mesonic case consists of the pseudoscalar meson H(v) and the 
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corresponding vector meson H*(v,e), since a spin rotation yields 

e W (iS v (e)^\H(v)) = (-i)\H*(v,e)) (14) 

where we have chosen an arbitrary phase to be (— i). The spin-symmetry 
doublets for baryons have been considered in and the general case, also 
valid for excited states, has been studied in . 

In the heavy-mass limit the spin symmetry partners have to be degenerate 
and their splitting has to scale as 1/mq. From the Lagrangian given above, 
one derives the mass relation for the heavy ground-state mesons up to terms 
of order 1/mg 

M H = m Q + A + (Ai + d H C 3 X 2 ) (15) 

where dn = 3 for the 0~ and dji — — 1 for the 1~ meson and C3 has been 
given in (|5p. Furthermore, the parameters A, Ai and A2 are given by 

(0\qj 5 h v \H{v))A = (0\q ivD j 5 h v \H{v)} (16) 
2M H \ l = (H(v)\h v (tD) 2 h v \H(v)) (17) 
6M H t\ 2 = (H(v)\h v (T lw iiy i iD v h v \H(v)), (18) 

where the normalization of the states is chosen to be (H(v)\h v h v \H(v)) = 
2Mh- These parameters may be interpreted as the binding energy of the 
heavy meson in the infinite mass limit (A), the expectation value of the kinetic 
energy of the heavy quark (Ai) and its energy due to the chromomagnetic 
moment of the heavy quark (A2) inside the heavy meson. 

All the parameters appearing in the mass relation are subject to renor- 
malization or suffer from ambiguities from renormalons, the latter subject 
is discussed in [T7]]. Hence quoting values for these parameters requires a 
procedure to be defined to deal with the ambiguities. 

The only parameter which is easy to access is A2, since it is related to 
the mass splitting between H(v) and H*(v, e). From the .B-meson system we 
obtain 

A 2 (w 6 ) = 1 -{M H * - M B ) = 0.12 GeV 2 , (19) 

and using the scaling @ we obtain the same value as from the correspond- 
ing mass splitting in the charm system. This shows that indeed the spin- 
symmetry partners are degenerate in the infinite mass limit and the splitting 
between them scales as 1/toq. 
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The other parameters appearing in ( |i5| ) are not simply related to the 
hadron spectrum. Using the pole mass for itiq in (fHf ), QCD sum rules yield 
for a value of A = 570 ± 70 MeV 0. More problematic is the parameter Ai; 
from its definition one is led to assume Ai < 0; a more restrictive inequality 

- Ai > 3A 2 (20) 

has been derived in a quantum mechanical framework in |T3] and using heavy- 



flavour sum rules 191. Furthermore, there exists also a sum rule estimate 



[ p0|] for this parameter: 

Ai = -0.52 ±0.12 GeV 2 , (21) 

which, however, has been critizised. A more extensive discussion of this issue 
is given in [f21"H . 

In the infinite mass limit the symmetries imply relations between matrix 
elements involving heavy quarks. For a transition between heavy ground- 
state mesons H (either pseudoscalar or vector) with heavy flavour / (/') 
moving with velocities v (v'), one obtains in the heavy-quark limit 

{H^'\v')\h { S ) Vh^\H^\v)) = f(w') Tr \u(^)TH(v)} , (22) 

where V is some arbitrary Dirac matrix and H(v) are the representation 
matrices for the spin structure of the heavy mesons 

JM^ (75^+0- meson 
H(v) = ^—^ UP+ 1- meson (23) 
[ with polarization e. 

The single form factor for these transitions, the Isgur-Wise function 
contains all the non-perturbative information for the heavy-to-heavy decay. 
Furthermore, heavy-quark symmetry fixes the value of £ at the point v — v' 
to be 

t{w' = 1) = 1, (24) 

since the current h\f'Th\/' is one of the generators of heavy-flavour sym- 
metry. The generalization of (22) to baryons may be found in [15] and to 
excited states in [ |i~6|j . 

The symmetries also place some restrictions on the corrections which 
may appear. In general, if explicit symmetry breaking is present those form 
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factors, which are normalized due to the symmetry, only receive second-order 
symmetry-breaking corrections. This general statement is the Ademollo- 
Gatto theorem E3|, which has been specialized to the case of heavy-quark 



symmetries by Luke [23 



2.2 Strategy of a HQET Calculation 

The relations ( |22|) and (Pl[) hold in the heavy-quark limit, and the machinery 
of HQET allows us to calculate corrections to (p2[) and (|24]) . 

In general there are two types of corrections. The short-distance correc- 
tions may be calculated in perturbation theory, based on the leading order 
of the 1/rriQ expansion of the Lagrangian. The logarithmic ultraviolet diver- 
gences in the effective theory correspond to logarithmic dependences on the 
heavy-quark mass tuq in the full theory, and renormalization group methods 
may be employed to perform resummations of these logarithms. In fact, the 
leading logarithmic corrections to bilinear currents are independent of the 
spin structure of the current. 

The second type of corrections are the power corrections of order l/mq, 
which in general involve long-distance physics and hence may in general not 
be calculated, but have to be parametrized. As an example, consider a matrix 
element of a current qTQ mediating a transition between a heavy meson and 
some arbitrary state \A). Using the expansion of the full QCD field (jl]), (0) 
and the corresponding expansion of the Lagrangian ([]), one has, up to order 

(A\qTQ\M(v)} = (A\qTh v \H(v)} (25) 

+^—(A\qTPADh v \H(v)) - i [ d^xiAYTiLAx^qTKMHiv)) + C(l/m 2 ) 
2m Q J 

where L\ are the first-order corrections to the Lagrangian as given in (|5|). 
Furthermore, \M(v)) is the state of the heavy meson in full QCD, including 
all its mass dependence, while \H (v )) is the corresponding state in the infinite 
mass limit. 

Expression (|25|) displays the generic structure of the higher-order cor- 
rections as they appear in any HQET calculation. There will be local con- 
tributions coming from the expansion of the full QCD field; these may be 
interpreted as the corrections to the currents. The non-local contributions, 
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i.e. the time-ordered products, are the corresponding corrections to the states 
and thus in the r.h.s. of (|25|) only the states of the infinite-mass limit appear. 



2.3 An Application: B — > D&£i/ 

As an application, we shall consider the weak transition B — ► D^*'£u, which 
is the most obvious, since both the b and the c quark may be considered as 
heavy. 

In general, the left-handed B — > D^> transition matrix elements are given 
in terms of six form factors 

(D(v')\c^b\B(v)} = [£+(</) K + <) + UvM - <)] (26) 

(D*(v',e)\c llM b\B(v)} = iy/m B m D *iv{y)e» aPp e* a v' p v p 
(D*(v',e)\c'j^ 5 b\B( y v)) = ^m B m D * 

Ui(y)(w' + l)e* - iA2(y)(e*v)v ll - Uiiy^v)^ 

where we have defined y = vv' . In the heavy-quark limit, these form factors 
are related to the Isgur-Wise function £ by 

&(2/)=£(j/)fori = +,V;Al,A3, 
&(j/) = 0fori = - A2. 



The normalization statement ( f24|) may be used to perform a model- 
independent determination of V c t> from semileptonic heavy-to-heavy decays 
by extrapolating the lepton spectrum to the kinematic endpoint v = v'. 
Using the mode B — > D^*'£u one obtains the relation 

1 dT G 2 

lim . — — = -^\V cb \ 2 {m B -m D *fm^\Ui{l)\ 2 . (27) 



47T 3 



In the heavy-quark limit the form factor £41 reduces to the Isgur-Wise func- 
tion and is unity at the non- recoil point; aside from \V c b\ everything in the 
r.h.s. is known. 

The corrections to this relation have been calculated along the lines out- 
lined above in leading and subleading order. A complete discussion may be 
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found in the review article by Neubert ||, including reference to the original 
papers. Here we only state the final result 



6/25 



1 + 1.561 



a s {m c ) - a s {m b ) 8a s (m c ) 



71 



3tt 



(28) 



154 27 18 25 J 



a, m z 



7T 1 — 2; 



In 2 



where we use the abbreviations 



Tfl r 



^(toO _ 

X . . , z 

a s (m b ) m b 

and to is a scale somewhere between m& and to c . 

The contributions in the square bracket originate from leading and sub- 
leading QCD radiative corrections. These include also the terms of order z n , 
which are short-distance contributions and hence may be calculated pertur- 
batively. 

The power corrections to the normalization are summarized in the cor- 
rection terms 5y m 2. The form factor £ai is protected by Lukes theorem, 
i.e. it does not receive corrections of the order 1/toq. Thus the first non- 
vanishing recoil corrections are of order (A/m c ) 2 , (A/m b ) 2 and A 2 /(TOfeTO c ). 
These contributions may only be estimated, since they need an input be- 
yond heavy-quark effective theory. There are various estimates for these 
corrections P^|-[P6] , which are compatible with one another; a very recent 



compilation of the various results yields 



■(5.5 ±2.5)%. 



(29) 



However, this is an estimate based on various assumptions; in fact the esti- 
mate of Si/ m 2 will be the final limitation for a model- independent extraction 
of V cb from exclusive decays. 



Adding all the corrections to the normalization, the value quoted in J27 
for the normalization is 



£ai(1) = 0.93 ±0.03, 



(30) 
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where the error of 6% is due to the uncertainty of the 1/m 2 corrections 
and the next-to-next-to-leading-order short- distance contibutions. This leads 
finally to a value of V c b, taking into account the latest data one finds P7 



\Veb\ = 0.040 ± 0.003. (31) 



3 The Heavy— Mass Limit for Inclusive De- 
cays 

Another important development in heavy-flavour physics was the formula- 
tion of the heavy-mass expansion for inclusive decays P- ||10|| , including even 



non-leptonic processes fl32fl . The main idea is to apply the operator-product 
expansion, making use of the fact that the heavy quark mass sets a large 
scale. This expansion involves operators with increasing dimension, the co- 
efficients of which are proportional to the appropriate power of 1/itlq. The 
mass dependence of the matrix elements of these operators may as well be 
expanded in powers of 1/mq using the machinery of HQET, and hence one 
may set up a 1/ mq expansion for inclusive rates and also for differential dis- 
tributions; generically the leading term of this expansion is the decay of a 
free quark. 

Applying this idea to the energy spectra of the charged lepton in inclusive 
semileptonic decays of heavy mesons, the relevant expansion parameter is not 
l/rriQ, but rather l/(mg- 2Ee); the denominator is thus the energy release of 
the decay. In almost all phase space the energy release is of the order of the 
heavy mass; it is only in the endpoint region that it becomes small and hence 
the expansion breaks down. This problem may be fixed by a resummation 
of terms in the operator product expansion, which strongly resembles the 
summation corresponding to leading twist in deep inelastic scattering. Anal- 
ogously to the parton-distribution function, a universal function appears, 
which determines all inclusive heavy-to-light decays. 

3.1 Operator Product Expansion 

The general effective Hamiltonian for a decay of a heavy (down-type) quark 
is given by 

H eff = QR (32) 
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where the operator R describes the decay products, e.g. 

R sl = ^V Qq Ta4 (1 - l5 )q (PeJ^l - J 5 )t) (33) 

for semileptonic decays. 

The inclusive decay rate for a heavy hadron containing the quark Q is 
then given by 



A' 



T oc Y,^)^\PB-q-Px)\(X s \H eff \B{v))\ 2 (34) 

' d A x(B(v)\H eff (x)Hl ff (Q)\B(v)) 

= 2 1m J d 4 x{B(v)\T{H eff (x)Ht ff (0)}\B(v)). 

The matrix element appearing in ([34]) contains a large scale, namely the mass 
of the heavy quark. The first step towards a 1 / ttiq expansion is to make this 
large scale explicit. This may be done by a phase redefinition as in (jl]). This 
leads to 

r oc 2 Im / d 4 xe- im Q vx (B(v)\T{H ef f(x)Hl ff (0)}\B(v)) (35) 

where 

7~(-eff — Q V R (36) 

with Q v from (|l|). In this way it becomes clear that a short- distance expansion 
is possible, if the mass tuq is large. The second step is thus to perform an 
operator-product expansion, which has the general form 

d'xe^ (B(v)\T{H eff (xM ff (0)}\B(v)) 



OO / 1 \ « 

Eh— C n+ M{B(v)\O n+3 \B(v)) 



n=0 



where O n are operators of dimension n, with their matrix elements renor- 
malized at scale /i, and C n are the corresponding Wilson coefficients. 

In a third step one removes the mass dependences from the matrix ele- 
ments by expanding the heavy quark fields appearing in the operators O n 
using (|l|) and (||), as well as the states by including the corrections to the 
Lagrangian as time-ordered products. 
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The lowest-order term of the operator product expansion is a scalar 
dimension-3 operator and hence it is either Q V ^Q V or Q V Q V . The first one is 
the Q- number current that is normalized even in full QCD, while the second 
may be related to the first via 



{iDf - (ivDf + -o„ v G» v 



^ + 0(1/7^). 

(37) 

where G^ u is the gluon field strength. Evaluating its contribution yields the 
free quark decay rate. 

All dimension-4 operators are proportional to the equation of motion 

04 oc Q v (ivD)Q v , and the first non-trivial contribution comes from dimension- 

5 operators and are of order of l/m%. For mesonic decays there are only the 
two parameters Ai and A2 given in (|i"7|) and fllSj) , which parametrize the 
non-perturbative input in the order l/trig. 



3.2 Inclusive Semileptonic Decays 

Applying the method outlined in the last paragraph to semileptonic decays, 
one finds for the decay B — > X c lv 



T(B^X c tv) = T b \V cb \ 2 x 
where 



1 + 



2w? I Vr/ifc 



1927T 3 



9A 2 . ( m c 



2m? ' 



and the two fj are phase-space functions given by 
fi(x) = 1 



x 2 + 8x 6 - x s - 2Ax 4 log x, 



f 2 (x) = 1 - -x 2 - 8x 4 + 8x 6 + -x 8 + 8x 4 logx. 
3 3 

From this one may read off the result for B — > X u lv^ 

Ai - 9A 2 



T(B -> XJu) = r b \v t 



ub 



1 + 



1m\ 



(38) 
(39) 
(40) 



(41) 



Expressions (|38D and fl4lD contain the leading non-perturbative corrections, 
parametrized by Ai and A 2 . However, before this may be confronted with 
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data, one has to apply as well QCD radiative corrections, which have been 
studied in detail [p8 



35 



The method of the operator-product expansion may also be used to obtain 
the non-perturbative corrections to the charged lepton energy spectrum. In 
this case the procedure outlined in the last paragraph is applied not to the 
full effective Hamiltonian, but rather only to the hadronic currents. The rate 
is written as a product of the hadronic and leptonic tensor 



dT 



G% 



Am 



B 



\V Qq \ 2 W^d{PS), 



(42) 



where d(PS) is the phase-space differential. The short-distance expansion is 
then performed for the two currents appearing in the hadronic tensor. Re- 
defining the heavy-quark fields as in (P and (|3|) one finds that the momentum 
transfer variable relevant for the short-distance expansion is mqv — q, where 
q is the momentum transfer to the leptons. 

The structure of the expansion for the spectrum is identical to the one 
of the total rate. The contribution of the dimension-3 operators yields the 
free-quark decay spectrum, there are no contributions from dimension-4 op- 
erators, and the 1/mf corrections are parametrized in terms of Ai and A2. 
Calculating the spectrum for B — > X c iv yields relatively complicated expres- 
sions, which may be found in |§-||. However, for the decay B — > X u lv the 
expression simpifies and is given by 



-- = 2, (3-2y) + -3-^ + 2y(6 + 5y)^ 



Ai + 33A 2 
3m 2 h 



tf(i-v) 



Ai 

3m? 



S'(l-y), 



(43) 



where y = 2Ei/m b is the rescaled energy of the charged lepton. 

Figure [I] shows the distributions for inclusive semileptonic decays of B 
mesons. The spectrum close to the endpoint, where the lepton energy be- 
comes maximal, exhibits a sharp spike as y —>■ y ma x- Close to the endpoint 
we have 

1,/J 26(1 -y-p) (44) 



T b dy 



Ai 



P 



(m Q (l-y)y 



3-4 
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Figure 1: The electron spectrum for free quark b —>■ c decay (dashed line), 
free quark b — > u decay (grey line), and B —>■ X c ev e decay including 1/mf 
corrections (solid line) with Ai = —0.5 GeV 2 and A2 = 0.12 GeV 2 . The figure 
is from ||. 
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where p = rn^jm^. This expression behaves like 5- functions and its deriva- 



tives as p — > 0, which can be seen in (43). This behaviour indicates a break- 



down of the operator product expansion close to the endpoint, since for the 
spectra the expansion parameter is not 1/rriQ, but rather l/(rriQ — qv), which 
becomes l/(mg[l — y]) after the integration over the neutrino momentum. 
In order to obtain a description of the endpoint region, one has to perform 
some resummation of the operator product expansion. 



3.3 The Endpoint Region E e ~ E s 



max 



Very close to the endpoint of the inclusive semileptonic decay spectra only 
a few resonances contribute. In this resonance region one may not expect to 
have a good description of the spectrum using an approach based on parton- 
hadron duality; here a sum over a few resonances will be appropriate. 

In the variable y the size of this region is however of the order of (A/mg) 2 
and thus small. In a larger region of the order A/mg, which we shall call the 
endpoint region, many resonances contribute and one may hope to describe 
the spectrum in this region using parton-hadron duality. 

It has been argued in p9[ that the 5-function-like singularities appearing 



in (|3) may be reinterpreted as the expansion of a non-perturbative function 



describing the spectrum in the endpoint region. Keeping only the singular 



terms of ( f43|) we write 

^=2y\3-2y)S(y), (45) 
T b dy 

where 

oo 

S(y) = Q(l-y) + Y / aJ^(l-y) (46) 

n=0 

is a non-perturbative function given in terms of the moments of the spec- 
trum, taken over the endpoint region. These moments themselves have an 
expansion in 1/mg such that a n ~ l/mg +1 , and we shall consider only the 
leading term in the expansion of the moments, corresponding to the most 
singular contribution to the endpoint region. 

Comparing (f43"D with (|45|) and ([46]) one obtains that 

a = Jdy(S(y)-e(l-y)) = (47) 
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ai = j y(S(y) - 9(1 - y)) = (48) 



where the integral extends over the endpoint region. 

The non-perturbative function implements a resummation of the most 
singular terms contributiong to the endpoint and, in the language of deep 
inelastic scattering, corresponds to the leading twist contribution. This re- 
summation has been studied in QCD [|30], [H| and the function S(y) may be 



related to the distribution of the light cone component of the heavy quark 
residual momentum inside the heavy meson. The latter is a fundamental 
function for inclusive heavy-to-light transitions, which has been defined in 


/(*+) = « ( B (v)\ KS(k+ - iD + ) h v \B(v)), (49) 

IrriB 

where k + = fc + ^3 is the positive light cone component of the residual 
momentum k. The relation between the two functions S and / is given by 

1 ? 

S(y) = — / dk+f(k + ) (50) 

TTIq J 

-mg(l-j/) 

from which we infer that the n th moment of the endpoint region is given in 
terms of the matrix element {B(v)\h v (iD + ) n h v \B(v )). 

The function / is a universal distribution function, which appears in 
all heavy-to-light inclusive decays; another example is the decay B — > X s ^y 
[ j3l"l , |18| , where this function determines the photon-energy spectrum in a 
region of order 1 / uiq around the K* peak. 

Some of the properties of / are known. Its support is — oo < k + < A, 
it is normalized to unity, and its first moment vanishes. Its second moment 
is given by a±, and its third moment has been estimated |18 , 25| . A one- 
parameter model for / has been suggested in [P0| , which incorporates the 
known features of / 

39 r 4 1 

/(M = -5T C 1 - x ? exp - - (1 - xf 6(1 - x), (51) 

where x = k + /A, and the choice A = 570 MeV yields reasonable values for 
the moments. In fig. |2| we show the spectrum for B — > X u lvi using the ansatz 
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Figure 2: Charged-lepton spectrum in B — ► X u iv decays. The solid line 
is (pE5D with the ansatz ([)!]), the dashed line shows the prediction of the 



free-quark decay model. The figure is from [30 



Including the non-perturbative effects yields a reasonably behaved spec- 
trum in the endpoint region and the 5-function-like singularities have disap- 
peared. Furthermore, the spectrum now extends beyond the parton model 
endpoint; it is shifted from E™'™ = m^jl to the physical endpoint EJ 1 ^ = 
Mb/2, since / is non- vanishing for positive values of k + < A = Mb — m^. 



3.4 Inclusive Non-leptonic and Rare Decays 

The same method as described above for the semileptonic decays has been 



ITCH and also to non-leptonic decays, see |32] for a recent 



applied to rare 

review. For simplicity we shall restrict the discussion here to a few simple 
examples. 

The effective Hamiltonian for the radiative rare decays B — > Xgj is given 
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by (B3) with 



R 



v/2 



327T 2 



/IV 1 



(52) 



where CV(/i) is a coefficient calculated e.g. in |33 



Going through the steps outlined above, one obtains for the inclusive rate 
of B — > Xg'-f, including the first non-perturbative correction [ lCf 



<xG\ 
16V 



ml\V ts V t * d \ 2 \C 7 (m b )\< 



1 + 



2m\ 



(Ai-9A 2 ) 



(53) 



Along the same lines one may study the inclusive decays B — > X s £ + £~ , the 
results for the lepton spectra and the total rates may be found in [|I(J . 

Similarly one may consider the total rates for inclusive non-leptonic de- 
cays. Neglecting for simplicity penguin contributions and CKM-suppressed 
decay modes, the effective Hamiltonian is given by 



Hp. 



Gi 



leff ~ 71 

with the two operators 



{c^o^ + cmom) 



01 = {&7^(1 - 7s)c}{u7 M (l - 75)4, 

2 = {^(1 - 75)4{"«7m(1 ~ 7s)c}, 



(54) 



(55) 
(56) 



where the braces denote the coupling to colour singlet, and the Cj are the 
corresponding Wilson coefficients, which may be found e.g. in ||34|1 . 
From this one obtains for the inclusive width for B 



B^X C 



G\m\ 
64vr 3 



■.nib 



l + -^.-AA,) 



- iSA 2 f 3 — 



1 



m b J 2m b 



. (57) 

where Ai/ 2 is a combination of Wilson coefficients appearing in the effective 
Hamiltonian 



At = Cl{m b ) + C 2 2 {m b ) + -Ci(m 6 )C 2 (m 6 ) 



(58) 



A, 



'-Ci(m b )C 2 {m b ), 
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and fi and fs are phase-space factors; fx is defined in (fTO|), while is 

/ 3 (x) = (l-x 2 ) 3 . (59) 

Equation (|57D is a QCD-based calculation of the inclusive non-leptonic 
width, and together with the corresponding expression for the semileptonic 
width (|38f) it gives us the lifetime of bottom hadrons and their semileptonic 
branching fraction, including the non-perturbative corrections in a 1/rriQ 
expansion. However, before one may compare the results with data, one has 
to take into account perturbative QCD corrections as well. These corrections 
have been calculated and are presented in a contribution to this conference 



4 Concluding Remarks 

The expansion of QCD in inverse powers of the heavy-quark mass has put 
heavy-quark physics on a model-independent basis. In particular, the symme- 
tries present in the heavy quark limit allow a variety of model-independent 
predictions for weak decay matrix elements. For mesons, the Isgur-Wise 
function ( ^2|) is the only non-perturbative input in the heavy mass limit. 

However, the corrections of order A/itlq introduce in general new form 
factors, i.e. an additional non-perturbative input is needed. Still a few rela- 
tions, like the normalization of certain form factors, are do not receive linear 
corrections, and the first subleading contribution is of order (A/mg) 2 . For 
tuq = m c = 1.5 GeV and for A = 500 MeV this gives a typical size of 
corrections in the ballpark of 10%, which is what is found, for example for 
the normalization of the form factor relevant for the V c b determination. In 
general, these corrections will be the final limitation for model-independent 
statements from HQET. 

These remarks apply in particular to the determination of from the 
exclusive channel B — > D*lv^ where the theoretical errors quoted above is 
about 6% and dominated by the uncertainties of the estimates of the l/^ig 
contributions, which need model input. This has to be compared with the 
determination of V c b from inclusive decays, which has been discussed in detail 
The inclusive width in the framework discussed above depends on 
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the quark masses, and superficially one finds a tjiq dependence. This would 
mean that even small uncertainties in the heavy-quark mass would have a 
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large effect on a V c b determination based on the inclusive width. However, it 
has been argued in [ 26 1 that the inclusive decays B — > X c lv will receive their 
major contribution from the kinematic region close to the non-recoil point; 
in this region the inclusive width depends almost linearly only on the mass 
difference — m c with only a weak dependence on + m c . The quark- 
mass difference is much better known than the individual masses; using the 
mass formula ( |15| ) only Ai enters. Based on these observations it has been 
argued in |26[ that a determination of V c b from inclusive decays may have 
a theoretical error of as low as 5%; the uncertainty here enters through the 
parameter Ai, which is at present only poorly known, but may be measured 



in the future from the inclusive semileptonic decay spectra [|37|]- Given the 



present situation both methods have comparable theoretical uncertainties 
and it remains to be seen for which of the two methods the uncertainties 
appearing at order will be better under control. 

Although there has been some theoretical progress in setting up a QCD- 
based calculation for inclusive widths, non-leptonic decays still remain a prob- 
lem. It has been noticed soon after the formulation of the I/tuq expansion 
for inclusive non-leptonic processes that the non-perturbative effects calcu- 
lated in this way are small, too small to explain the experimental data on the 
inclusive semileptonic branching fraction of B mesons. However, there are 
perturbative corrections as well, which have been calculated recently, taking 
into account a non-zero mass for the quarks in the final state |35], [3(J . These 
corrections are substantial only in the channel b — > ccs and hence yield an 
enhancement charm production in B decays that is not supported by present 
data. Thus the problem of the semileptonic branching fraction still persists. 

The difficulty seems to be the calculation of the inclusive non-leptonic 
width, and not the semileptonic one. This is supported by another problem, 
which is the lifetime of the baryon. Based on the 1/tuq expansion one 
would conclude that the lifetime should be slighly smaller than the B 
meson lifetime, t a ~ 0.9r B [O. This is not supported by recent data, 



indicating that r\ b ~ 0.7tb where the experimental error is 15% |38 |. The 
situation in the charm system is even worse, here the lifetime differences 
are substantial, t\ c ~ 0.5r D o and ta c ~ 0.2r D ±. This indicates that the 
1 / rug expansion for inclusive non-leptonic decays is not yet understood and 
the problems have been recently summarized in [RjJ. Unlike exclusive non- 
leptonic decays, which still may be described only in a model framework, the 
description of inclusive non-leptonic decays is based on QCD and the above 
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problems certainly deserve further study. 
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